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Abstract 

For decoherence processes induced by weak interactions with the 
environment, a general quantum channel with one noise parameter 
has been formulated. This channel is called low-noise channel and 
very useful for investigating the parameter estimation in the leading 
order. In this paper, we formulate the low-noise channel with multi- 
ple unknown parameters in order to address the simultaneous achiev- 
ability of the Cramer-Rao bound for the parameters estimation. In 
general, the simultaneous achievement of the Cramer-Rao bound for 
multi-parameter estimations suffers from non-commutativity of opti- 
mal measurements for respective parameters. However, with certain 
exceptions, we show that the Cramer-Rao bound for output states of 
dissipative low-noise channels can be always attained in the first order 
of the parameters as long as D < — 1, where D and N denote the 
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number of the parameters and the dimension of the system, respec- 
tively. This condition is replaced hy D < N'^ — 1 if it is allowed to 
set the entanglement with ancilla systems in its input state and to 
perform the non-Iocal measurement on the composite system. 
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1 Introduction 



Decoherence processes driven by weak interactions with the environment 
can be described in a unified way by the low-noise channel formulation pro- 
posed in Ref . [1] , which has a wide range of applications for the estimation 
of small noise parameters on quantum channels. An important example of 
its application is the noise estimation in quantum computing [2] in which 
small noises induced by the environment can lead to a serious obstacle to 
its implementation. One of the best way to eradicate the noise effect is to 
apply an appropriate error correcting code to the system. In this scheme, one 
first needs to know properties of the noises to design the code. However, be- 
cause the noises are generated by very complicated many-body interactions, 
it is difficult to calculate the properties theoretically from first principles. 
Therefore, it is effective for finding out the properties to measure the noise 
experimentally and to estimate them from the experimental data. With this 
method, the estimation of such very weak noises may suffer from large am- 
biguity of experimental errors. Thus the estimation theory of the low-noise 
channel will play a significant role in order to improve the performance of the 
estimation. Another example to which the estimation theory of the low-noise 
channel can be applied is found in rare processes of the elementary particle 
physics [3] . Many theories of new particles has predicted the existence of very 
weak interactions which generate rare decays and rare reactions of particles. 
However, in real experiments, the number of signals indicating new physics is 
generally small, compared with that of standard-model signals. Thus, siding 
with the standard model, the evidential signals can be regarded as low-noise 
background datum. The estimation theory of the low-noise channel is ex- 
pected to improve accuracy of estimation of such tiny reaction rates and to 
assist in the finding of new physics. 

An estimation theory for the low-noise channel has been formulated for 
one noise parameter, termed a low-noise parameter, which is assumed to be 
unknown and very small [1]. The increase of the Fisher information due to 
the prior entanglement with non-local output measurement, which is called 
the ancilla-assisted enhancement, has been analyzed for the ancilla-extended 
version of the low-noise channel [1]. The low-noise channel has been also 
extended to the channel on the many-body system [4], and this work shows 
that the maximum of the Fisher information can be attained by a factorized 
input state. 



1 



In this paper, we address the estimation problem of multiple low-noise 
parameters e = (e^, ■ ■ ■ , e^), where D is the number of the parameters. We 
assume that all of are non-negative and take very small values of the 
same order. In later analysis, 0(e) denotes the order of e^, and the orders 
of products of the parameters, 11^=1 (^^)"'') denoted by O (e^i'*' '+"^). 
Due to the multi-parameter degrees of freedom, a large variety of physical 
phenomena can be analyzed by use of this low-noise channel. 

The low-noise channel are dependent on the unknown parameters e which 
are estimated by measurements. Consider an input state p^n, which is inde- 
pendent of e, into the low- noise channel F^. The output state p = Pout{^) is 
obtained by 

Po«t(e) = rjpi„]. (1) 

If the purity of the output states is decreased, then the low-noise channel 
is called dissipative. In what follows, we focus on the dissipative low-noise 
channel. 

To estimate the low-noise parameters, a POVM measurement with D 
measurement values are performed on the output state Pouti^)- The mea- 
surement values are regarded as estimation values for and this POVM 
is called an estimator. When the expectation value of the estimation value 
is equal to the true value of the parameter for any p, the estimator is 
called locally unbiased. In the present paper the locally unbiased estima- 
tor is defined in the vicinity of e = 0. Note that by virtue of the Naimark 
theorem [5, 6], the POVM measurement is equivalent to a measurement of 
commutative observables (projection- valued measures) in a composite system 
comprising the original system and an ancilla system. 

One of the fundamental quantities in the low-noise estimation theory is 
the Fisher information matrix of an output state of the low-noise channel 
Ff. It is well known as the Cramer-Rao inequahty that the inverse of the 
Fisher information matrix is a lower bound of the mean-square-error matrix 
of the locally unbiased estimators [5, 6]. In this paper, we focus our attention 
on the achievement problem of the lower bound for the dissipative low-noise 
channel. In the case with one low-noise parameter, the bound can be attained 
by using certain estimators. In contrast, the bound cannot be always achieved 
if D > 1. Even if a certain unbiased estimator is optimal for the parameter 
estimation of e^, the estimator is not always optimal for other parameters 
e^, • • • ,e^, because of the non-commutativity of optimal estimators. It is 
an important open problem to reveal the necessary and sufficient condition 
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for the attainability of the bound for the muhi-parameter estimation. The 
solution for general channels has not been found to date. Recent progress 
has been reported on the multi-parameter estimation for special families of 
quantum channels, in particular, SU(N) channels [7, 8, 9] and a generalized 
Pauh channel [10]. 

Our main results are related to the attainability problem for the esti- 
mation in the dissipative low-noise channel. To explain the results, let us 
diagonalize the output state p{e) as 

N-1 

p(e) ^Po(6)|0(6))(0(6)| + J2pn{e)\n{e)){n{e)\, 

n=l 

where the first eigenvalue Poie) is 0(e°), the other eigenvalues Pn^oi^) are 
0(e), and denotes the dimension of the Hilbert space of the quantum 
system (the reason why this expression is available is shown in Sec. 4). The 
Fisher information of the output state for the channel has a leading term of 
order of 0(e~^). We first prove in the case with D < N — 1 that all output 
states of the dissipative low-noise channel have locally unbiased estimators 
by which the lower bound is attained in 0(e) assuming that 



det 



N-1 



,n=0 



^0 (2) 



in the vicinity of e = 0, where stand for d/de^. In order to explain 
the meaning of Eq. (2), let us consider a map {\fpo{e), • • • , ^JpM-l{,^)) which 
generates a manifold by running the D parameters of e. Clearly, ilD > N — 1, 
Eq. (2) does not hold because the target space defined by 

N-1 

(Vpo, ■■■ , VPn-i) Pi ^ 0, = 1 

i=Q 

has N — 1 dimensions and the parameterization of the manifold by e becomes 
degenerate. On the other hand, when D < N—1, it is noticed that imposition 
of Eq. (2) is not so difficult. If the dimension of the manifold is equal to D, the 
relation in Eq. (2) automatically holds. Therefore, the condition in Eq. (2) 
is not anomalous and can be satisfied for many standard low-noise channels. 

Next, we point out that, although the lower bound is not always achieved 
even in 0{e^) when D > N — 1, the entanglement with ancilla systems can 
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make the lower bound attainable again in 0{e) when — 1 > D. As is ex- 
pected from the fact that the dissipative low-noise channel comprise a rather 
general class of quantum channels, these results are very surprising because 
they reveals that an enormous number of channels can overcome possible 
non-commutativity among the optimal estimators to attain the Cramer-Rao 
lower bound. 



2 Brief Review of State-Parameter Estima- 
tion 



We start with a brief review of a parameter estimation theory of quantum 
states, since the channel estimation theory is based on the state estimation 
theory. Let us consider a state p dependent on D unknown parameters O^^, 

p^p{e)^p{e\---X). (3) 

The parameters 9^ are estimated by measuring the state p with a positive 
operator- valued measure (POVM) {n^|x = • • • ,x^)] which satisfies 

> 0, j XlJ^'x = 1, (4) 

by definition. The POVM Ti^ is called the estimator for the parameter esti- 
mation. When the condition 

j x'^Tr[n^p(^)]d^x = e^" (5) 

is satisfied, the estimator is called unbiased. For the unbiased estimator, the 
mean-square-error matrix V^^ is defined by 

y^^[n,] := j {x^ - e^){x'' - r)Tr[n,p(^)]rf^a;. (6) 

As is well known, there exists a universal lower bound of the error matrix 
V^^ , called the Cramer- Rao bound [5, 6], as 

V>J-\ (7) 
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where J stands for the Fisher information matrix for p defined by 

J^,= iTr[p(L^L, + L,L^)], (8) 
and L^'s are the symmetric logarithmic derivative (SLD) defined by 

Ll^W^ ^p^l{L^P + pL,). (9) 

Eq. (7) is called the Cramer-Rao inequality and implies that for an arbitrary 
£)-dimensional vector it, the following relation is always satisfied, 

'uVTi > itJ-^Tt. (10) 

Note that the equality is not always attainable. This aspect is related to 
the non-commutativity of quantum measurements. When D = 1, measure- 
ment of the SLD operator Li is optimal for 9^ estimation, that is, the equality 
V — is satisfied. When D — 2, each measurement of the SLD operator 

is optimal for 6'^ estimation independently However, the observables Li 
and L2 do not generally commute to each other. 

This non-commutativity of SLD's makes achievement conditions of the bound 
nontrivial. For the pure-state case {p{0) — \(j){9)) {(j){9)\) , it has been proven 
in Ref. [11] that the bound can be always attained. However, generally, the 
necessary and sufficient condition for the bound achievement has not been 
known. 

The state estimation theory explained thus far can be applied to the chan- 
nel estimation done as follows. Let us consider a ^'^-parameterized channel 
Tg. The channel generates an output state p{9) = Tg[pin] from a input state 
Pin independent of 9^^. The parameter on the channel is estimated by use 
of the state estimation theory for the output state p{9). 

The Fisher information matrix J[pin] for the output state changes as the 
input state pin changes. Hence, it is important to find input states which 
optimize J. For this problem, there is a useful theorem [12] related to choice 
of the input state. There exists a pure input state |0)(</'| for which the Fisher 
information matrix J is the maximum over all input states pin, that is, for 
an arbitrary vector u ,the following relation holds. 

lIJ[\cf>){cf>\]lt>lIJ[pin]lt. (11) 
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As was seen in Eq.(7), the mean-square-error matrix V is lower bounded by 
the inverse of J, therefore, without loss of generality, we may concentrate on 
pure input states in later discussions. 

As mentioned above, the Cramer-Rao bound for the output states cannot 
be always attained and the necessary and sufficient condition for its achieve- 
ment in channel estimation is a crucial open problem. For unitary channels, 
the achievability problems of the bound can be completely solved by use of 
the result in Ref. [11], since output states of the channels for pure input states 
remain pure, and thus the estimation problem of the channels is reduced to 
that of the output pure states. To date, the solution of the lower-bound 
achievement has not been found for general channels. 

3 Multi-Parameter Low-Noise Channel 



In this section, we define the low-noise channel Fg with unknown multi- 
parameter e. This channel is a natural extension of the low-noise channel 
with one parameter first introduced in Ref. [1]. The definition is given by a 
Kraus representation as 

Pie) = r,[p] = Yl Ba{e)pBi{e) + Y.''Y1 C,a{e)pCl{e), (12) 

a A'=l a=l 

where the Kraus operators satisfy the following four conditions. 

(i) The channel is a TPCP map: 

D 

E Bi{e)B^{e) + E E CUi^)C,a{e) = I5. (13) 

a n=l a=l 

(ii) Ba{e) is analytic at e = 0, giving the power series expansion 

B^{e) = K^ls- E N^:^^-^'^\e^f'---{e^r'^, (14) 

(ni,... nD)^(0,- ,0) 

in the neighborhood of e = 0, where and N^'^'"^'^^ are coefficients and 
operators, respectively, independent of e. 
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(iii) Ka satisfies 

a 

(iv) C^a(e) is analytic at e = 0, giving the power series expansion 

C,a{e) = + Yl M(:™)(e^)"i • • • (e^)"-, (15) 

(ni,... nD)^(0,-,0) 

where M^^ and M^m '"'^"^ arc operators independent of e. 

The condition (i) is a natural characteristic of physical channels, because 
all physical channels are TPCP maps and any TPCP map has Kraus repre- 
sentations. The conditions (ii) and (iv) simply imply that the channel shows 
nonsingular behavior near e = 0. Therefore, taking proper limits of weak 
coupling with the environment, physical processes induced by the environ- 
ment can always be described by this low-noise channel. From condition 
(iii), the channel automatically reduces to the identical channel in the limit 
of vanishing the parameters: 

lim Fg = ids-1 (16) 

e— »+0 

where ids stands for the identity channel on the system S. This shows that 
the parameters can represent noise parameters caused by the environment. 

The second term of the right-hand-side of Eq. (12) describes a dissipative 
effect by the environment in the lowest order of e and plays a crucial role in 
the later discussion. Without loss of generahty, we may assume that each of 
C^a is not proportional to others. If not, the second term can be rewritten in 
the form of Eq. (12) with a smaller number of independent C^a(e) by using 
coordinate transformation of e^. For the sake of convenience, let us denote 
K the total number of C^a(e) as 

D 

ij.=i 

In addition, we redefine the Kraus operators as 

N,a = N^r-'\ N,^ = 7V(oio-o)^ • • • , iVzPa = N^r-^^\ (17) 
Then Ba{e) is rewritten as 

D 

B^{e) ^K^l-Y, ^^N^a + 0{e'). (18) 
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The TPCP condition (Eq. (13)) in the lowest order of the parameter e is 
given by 

a a 

Due to this relation, the operator i<^a-^na can be broken down into a sum 
of two terms such that 

E <^Ma = ^ E ^l^^<^ + ^^M, (20) 
a a 

where arc Hermitian operators. Using Eqs. (12) and (20), the derivative 
of p(e) at e = is evaluated as 



9mP(0) = E 



1..., 1 



M^„p(0)Mt„ - -MiM^„p(0) - -p(0)MiM^„ 



a=l 

-i[H,,pm. (21) 

The first term of the right hand side indicates a noise effect with decoherence 
induced by the environment. On the other hand, the second term expresses a 
unitary time evolution of the system S. Therefore, if the first term vanishes, 
the purity of the output states does not change in 0(e). In the later dis- 
cussions, we assume that the operators M^a are non- vanishing and thus, in 
general, the second term does not disappear. This assumption implies that 
decoherence by the environment inevitably takes place in the system S. In 
this case, the low-noise channel is called dissipative. 



4 Channel Estimation 



We now discuss the attainabihty of the Cramer-Rao inequahty for the 

multi-parameter estimation in the low-noise channel. Let us consider a pure 
input state for the low-noise channel Fg. The output state p(e) = 

re[|0)(0|] can be diagonalized as 

N-l 

p{e) = Y.Pn{e)\n{e)){n{e)\ (22) 

n=0 
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with an orthonormal basis {|n(e))}. Because of a property of the low-noise 
channel (See Eq. (16)), wc arc able to impose the following boundary condi- 
tions, p(0) = 10) (01, po{0) = 1, and PnM^) = 0> which imply |0(0)) = 

The deviation operator Sp{e) = p(e) — |0)(0| between the input state and 
the output state is evaluated by use of Eq. (21) as 



<5p(e) = ^£'^a^p(0) + O(62) 

/J, a 

-iJ2^nH,: |0)(0|] + O(6^). (23) 

Here, let us define an — 1 dimensional matrix by 

A(e) = [(n(0)|5p(e)|n'(0))]. (24) 

Note that all of the orthonormal vectors |^(0)) (n = 1, ■ ■ ■ , — 1) are or- 
thogonal to 10) by definition. Then we see that the matrix A(e) is rewritten 
using Eq. (23) such that 



A(e) 



5]e'^(n(O)|M^„|0)(0|MtJn'(O)) 



fj,a 



+ Oie'). (25) 



The leading term of A(e) is 0(e), because, as mentioned before, the Kraus op- 
erators M^a do not vanish due to the assumption of the dissipative low-noise 
channel. Thus the matrix A(e) possesses 0(e) eigenvalues. Some eigenvalues 
of A(e) may be higher terms as O(e^) or exactly zero. However, the presence 
of such eigenvalues does not affect the following analysis at all. 

The matrix A(e) becomes a diagonal matrix with the basis {|^(0))} in 
0(e) and the — 1 eigenvalues arc 6pn defined by Spn{e) = Pn{^) — Pn(0), 
that is, 

Sp{e)\n{0)) = 6pn{e)\n{0)) + 0{e^), (n = 1, • • • , AT - 1). (26) 

It is worthwhile noting that the eigenvalues 5pn{e) are unchanged even if we 
use a different orthonormal basis |m(0)) := '^^Umn 1^(0)) in the definition 
of A(e), where Umn is a unitary matrix. Thus, when one calculate 5pn{e), the 
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orthonormal basis {|?T'(0))} can be arbitrarily fixed in Eq. (26). In addition, 
for < — 1, non- vanishing 5p,„(e)s can be obtained by solving another 
eigenvalue equation (See Appendix). Using the eigenvalues Spn{e), we obtain 
another expression of p(e) such that 

N-l 

p{e) = |0)(0|+5^e'^5] (a>„(e)|n(e))(n(6)| + \d,0{e)){<P\ + \<P){d,0{e)\)+O{e'), 

jj, n=0 

(27) 

where |9^0(e)) means (9^ 1 0(e)). From Eq. (27), the following relation is 
straightforwardly obtained, 

{n{e)\d^p{e)\m{e)) = d^5pn{e)5nm + 0(e), (28) 

for n, m 7^ 0. 

Now we describe the Fisher information matrix for the low noise channel. 
From Eq. (22), the SLD's are calculated [5] as 

{n{e)\L,{e)\m{e)) = I {n{e)\d,p{e)\m{e)) . (29) 

Using this, we have the Fisher information matrix J(e) — [J,j,v{€)] as 

J,A^) = ^(n(6)|a^p(6)|m(6)) I {m{eMp{e)\n{e)). (30) 

nm Pn{e)^Pm{e) 

If all of the eigenvalues 5p„(e) are 0(e), the Fisher information matrix J(e) 
is evaluated from Eq. (28) as 

J ^,^JfaM,^mpM^O(e»), (31) 

If some of the eigenvalues 6pn{^) are O(e^) or exactly zeros, the terms asso- 
ciated with the O(e^) eigenvalues in the sum of Eq. (31) can be neglected. 
Taking account of a fact that the leading order of J(e) is 0(e~^), let us 
define a divergent part of J(e) by 

N-l 



TdivfA ._ df,6pn{e)djpn{e) 



n=l 
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It is worthwhile noting that the classical Fisher information J^^{e) defined 

by 

W ■■= E ^^^^4^^ = 4X;9.V^5.^ (33) 



n=0 



coincides with J*^(e) in 0(e ) 



The classical information matrix J^^\ 
sional manifold M embedded in an - 

■ , ^/PN^) \Pi ^ 0, ^^0^ 



n=0 

, FO(e°). (34) 
e) is an induced metric of a D dimen- 



1 dimensional space P defined by 
Po, ■ ■ ■ , ^JVn-\) \Pi ^ U, l^i^Q-pi = l}. When (e^ ■■■ ,e^) locally 
has a one-to-one correspondence with a point on the manifold M, the matrix 
J^^(e) has its inverse matrix J'^'^'^{€) which is 0(e). If D > — 1, some pa- 
rameters of e are redundant and the parameterization of Pn{^) by e'* becomes 
degenerate, that is, 



det 



'N-l 



n=0 



0. 



Hence, the inverse matrix J'^'^^[e) does not exist. On the other hand, when 
D < N — 1, the inverse matrix J'^''^^{e) exists for Pn(e) non-degenerate in the 
neighborhood that e = 0. Then, because of this fact and Eq. (34), the inverse 
matrix J'^™-M'^(g) of J^^^(e) also exists and its order is 0(e) just as J'^'^^(e). 
In later analysis, we assume non- degeneracy of the parameterization; 

"7V-1 



det 



,n=0 



7^0 



(35) 



for D < — 1. In this case, the quantum Fisher information J(e) also has 
its inverse matrix J~^{e) — [J'^^{e)\, whose order is 0(e), and satisfies that 



J^"(e) 



(36) 



Now we shall prove that a locally unbiased estimator, which attains the 
Cramer- Rao bound in 0(e), can be explicitly constructed. Firstly, let us 
introduce commuting Hermitian operators A^{e) given by 



N-l 



n=l 



'9^'^Pn(e) 
(^Pn(e) 



|n(e))(n(e)|. 



(37) 
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If some of 5pn(e) are O(e^), the corresponding terms are neglected in the sum 
of Eq. (37). Their contravariant operators A^^{e) are also defined by 

A^ie) := J2 J'^'''-'""{e)A,{e). (38) 

Note that the order of yl^(e) is 0(e°). This operator can be written by the 
spectral decomposition as 

A'^{e) = J2<i^)Pn{e), (39) 

n 

where Pn{^) are projection operators satisfying Pn{^) = 1- Here we define 
a POVM operator by 

N-l D 
n=l /U=l 

We adopt as an estimator for a measurement to estimate e. Then it is 
proven from Eq. (28) that the estimator satisfies 

J x''T:t[U^p{e)]d^x = Tr[>l''(0)p(e)] 

= e^ + Oie^). (41) 

Thus the estimator Ila, is locally unbiased for the e estimation in the vicinity 
of e = 0. It is also shown that the estimator locally attains the Cramer-Rao 
bound as follows. Using Eq. (37) and Eq. (38), we have 

^Tr[p{e){A^{e), A-^ie)}] = J*-'^''(e) + 0{e'). (42) 
On the other hand, the mean- square-error matrix V'^'^ for is evaluated as 

V- = /(.--.^)(."-OTV[n.pWid-. 

= iTVW£){>l''(6), A-'{e)]] + 0(e'). 

Consequently, we see that the Cramer- Rao bound is actually satisfied in 0(e): 

= 7/^^(6) + 0(e2). (43) 
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This is a main result of this paper. This imphes that the lower bound can 
always be achieved in the estimation of the dissipativc low-noise channel 
for D < — 1 if it satisfies Eq. (35). Because, generally, multi-parameter 
estimation cannot attain the bound, this result is very significant. 

In the case with D > N — 1, we are able to give another important discus- 
sion. The inverse matrix J~^{e) of the quantum Fisher information matrix 
may exist in this case, even though the classical Fisher information matrix 
J^j,(e) does not have its inverse. However, even if the inverse matrix exists, it 
is not sufficiently suppressed in general. Its components are generally 0(e°), 
not 0(e). Moreover, the Cramer- Rao bound cannot be always attained even 
in 0(e°). The disadvantage can be remedied by using the ancilla-extension 
of the low-noise channel [4]. In fact, the bound becomes attainable again 
when D < N'^ — 1. This is the second result of this paper. 

To see this, let us consider an ancilla system A which has a Hilbert space 
with the same dimension N as that of the original system S. Suppose that 
the dissipative low-noise channel is extended to Fg^id^, where id a stands for 
the identity channel on the ancilla system. Entangled states of the composite 
system S + A are available as input states for the channels and a collective 
measurements is performed for the output states in order to estimate the 
D parameters e'^. In this case, the dimension Ns+a of the Hilbert space of 
5" -I- ^4 is larger than that of S because Ns+a — N"^- Therefore, when the 
number of the parameters is less than Ns+a {D < Ns+a — 1), the above 
analysis shows that the Cramer- Rao bound can be attained in 0(e) assuming 
that the ancilla-extended channels satisfy Eq. (35) for their output states. 
Consequently, when N < D < N"^ — 1, the ancilla- assisted enhancement 
effect is able to make the Cramer- Rao bound achievable in 0(e). 

5 Examples 

Let us consider a two-parameter dissipative low-noise channel Fg for a 
system S with N — 3 and D — 2 which is a simple example satisfying 
K < N — 1, and focus on the case with Ki — K2 — 1. Suppose that an 
input state |0) go through the channel. We apply the method described in 
Appendix to solve the eigenvalues of the Sp±(e) of the output state. In this 
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present system, the matrix A(e) defined by Eq. (55) can be written as 



A(6) 



where denotes M^i with /i — 1,2,3. The following eigenvalues is obtained 
by solving Eq. (54), 

5p±{e) 
_ 1 
" 2 

= 0(6), 



where (5M„y are variance- matrix elements defined by 



5M21 SM22 



The inverse of the Fisher information matrix can be explicitly calculated as 

(e^)^ 5Mu det SM + {e^f e^^Mss [3SM12SM21 - 26MnSM22] + {e'^f SM^^ 



J 
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{SMuSM22 - SM12SM21) {e^Mu - e^SM22) 
{e^f 5M22 det SM + {e^f e^SM^ [36M12SM21 - 2SM11SM22] + {e^f SMf^ 



{5MuSM22 - SM12SM21) (eMMn - e^t^Mss)^ 

5M12SM21 e^SMii + e^5M22 

5Mii5M22 - 5Mi25M2i (eMMn - e^^Maa)^' 
Hence, for input states which satisfy 



J 



12 



J 



21 



(44) 



5MuSM22 - SM12SM21 
n 



0(e°) 



with := e^/e, the inverse matrix behaves as 0(e). Thus, the Cramer- 
Rao lower bound can be attained in 0(e), since one can make an locally 
unbiased estimator to attain the bound by following the procedure in the 
previous section. 

Next, let us give another example for N = 2 and D = 2 to show the 
ancilla-assisted improvement as mentioned in Sec. 4. Now suppose that an 
input state pj„ goes through a dissipativc low-noise channel defined by 

Te [pin] = (1 - e^ - e^) Pin + e^CTxPin'^x + ^^'^zPin'^z, 
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where e'^ are unknown low-noise parameters and and ay are the Pauh 
matrices. Its output state p(e) can be expressed by a Bloch representation 

as 

p(e)=f,[p,„] = - + -y-a, (45) 

where y is a Bloch vector and a — (ctj;, Uy, cr^). The SLD for the output states 
is exactly solved and written as 

= lo,t, + ln-^, (46) 

where 



21-|y|2' 
and 

The Fisher information matrix J^^, is directly calculated as 



ld,\y\'dM 



2 



Defining 5p = 4(1 — |yp) which vanishes when e = 0, the inverse matrix J^*^ 
is calculated as 

tj — 



J 



\diW\d2y\'^ - {diy- ^2^)2 + ^ \d2\y\'^diy- di\y\^d2y\ 



2 ' 



\d,y\'\d2y\' - (9iy • d^y)' + ^ \^2\y\'^^y- dM^d^yl 



2 ' 



diy ■ d2y + j^di\y\'^d2\y\ 



2 



J,, ^ . . .-^ ^^g^ 

|<9i^P|(92^|2 - (^i^- d2yy + ^ \d2\y\^diy- di\y\'^d2y\ 

By solving the eigenvalue equation of J~^, we can show that one of the 
eigenvalues is 0(e), but the other is 0(e°), namely, 



diag[J ^] 



O(e0) 
0(e) 



(49) 
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Therefore, due to the 0{e'^) lower bound, it is impossible to estimate the 
parameters with 0(e^/^) precision. This property is briefly illustrated when 
e = at which the inverse matrix can be written as 



J-\0) 



1 



d2\y\'^ 
-di\y\ 



with $ := \d2\y\ diy — di\y\ d2y\ ■ It is easily seen that J has a zero 

rp 

eigenvalue with an eigenvector [ (9i|^p i92|yp ] , that is. 



J-\0) 





y\ 


_d2 


y\ 



0. 



It turns out in 0{e°) that the vector [ d2\y\'^ ]^ corresponds to the 

eigenvector with the eigenvalue of O (e) and the other eigenvector [ ^2 1 y| ^ — | j 
has the eigenvalue of 0(e°). 

In order to attain the Cramer- Rao bound of 0(e), let us consider an 
ancilla-extension of the channel, (g) id a, and its entangled input state l^') 
described by 



I*) 



-)), 



(50) 



where cr^li 



±|±). Then, the output state is given by (r^ id^)[|*)(*|]. 
In order to evaluate the eigenvalues Spn{e), let us define three vectors 





71^1 + ^^ 


|2) = 




|3) = 





(51) 



Note that these vectors are mutually orthogonal and all of them are orthog- 
onal to the input state |^). Thus it is possible to represent the matrix A(e) 
(See Eq. (24)) with a basis {|1), |2), |3)} as 



A(6) 



e2 
eV2 eV2 
eV2 eV2 



(52) 



16 



By estimating eigenvalues of A(e), it turns out that all the three eigenvalues 
behave as 0(e). In fact, the three eigenvalues are evaluated as 



5pi{e) = e^ 
5p2(e) = e\ 
Spsie) = 0. 



(53) 



Substituting 5p„(e) into Eq. (31) gives the Fisher information matrix such 
that 

i 



+ 0(e') 



Here the term d^6p-i{e)d^6ps{e)/Sp-i{e) in Eq. (31) has been set to zero. It 
can be verified that a mapping onto the manifold M defined by 

(VPo, VPi^ VP2, Vps) := (^Vl - <^Pi(e) - Sp2{e) - Sp3{e), y/Spi{e), y/Sp^, y/Sp^^ 
is non-degenerate; 

Yfn=o9lVPn{(^)^l^/pn{e) J2l=o di^ Pn{e) 82 y/ Pn{e) 



det 



^0. 



This non-dcgcncracy guarantees the existence of [J'^"]. The matrix [J'^"] is 
suppressed as 0(e). In fact, [J'^''] is given by 



[J' 




e 



-2 



Calculation of the projective operators {|n(0))(Ti(0)| \n — 0, 1, 2, 3} is simple 
and results in 



|0(0))(0(0)| = |*)(*|, 

|i(o))(i(o)| = i(l + +)-!--))((+ + 1 -(--I), 

|2(0))(2(0)| = ^(1 + -) + I -+))((+-! + (- + I), 
|3(0))(3(0)| = ^(1 + -)- I -+))((+-!-(- + I). 

Adopting the estimator, the ancilla-extension of the channel really makes the 
Cramer- Rao bound achievable in 0(e). 
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Appendix 

Here, we shall show that the non- vanishing (5p„(e) satisfy 

det[5p„(e)lK-A(e)] =0, (54) 

for K < N — 1, where the matrix A is a x Hermite matrix defined by 

A(e) = [A'^'''\e)] 

"^((01 {Ml - {mla\<f>)) (M., - (0|M.,|0)) 10)) V^] .(55) 

This means that one can evaluate 6pn{^) by solving Eq. (54). 

The proof of Eq. (54) is as follows. We begin with a function /(p) defined 

by 

f(p) := det [plAT-i - A(e)] . (56) 
This function can be expressed by products of Tr [A(e)*^] as 

f{p) = p^-i-p^-2TrA(e) + --- 

N-l 

= p^-^ + P'^"'""^n (Tr [A(e)"] , Tr [A(e)"-i] , ■ ■ ■ , Tr [A{e)]) . 

n=l 

(57) 

where {C„(-)} denote some coefficients expressed as a function of Tr[A(e)*^] 
with k — 1,2, - • ■ ,n. In Eq. (57), the trace of A(e)'^ is calculated by use of 
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the relation '^^=1 — Iat — \<f>){<P\ ^ 

Tr [A(e)T 

••■M^.,.|0)(0|M;,^,(liv-|0)(0|)] 
• • • V^{(l)\My,^,_, {In - |0)(0|) M^.„. 

= IV [A(6)T . 
Thus, for X < A?" — 1, we obtain 

N-l 

f{p) = p^-^ + 5^p^-^-"C4Tr[A(er],Tr[A(eri],---,Tr[A(6)]) 

n=l 
K 

ft=i 

= p^-^-^ p^ + f]p^-'^C„(lV[A(er],IV[A(6r-i],---,IV[A(6)]) 

n=l 

= p^-^-^detblx-A(6)]. 
Consequently, applying this result to the eigenvalue equation for A(e) gives 

det [6pni€)lK - He)] =0, 
which completes the proof. 
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